Planar array of semiconducting nanotubes in external electric field: Collective 

screening and polarizability 
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We study theoretically the charge separation in a planar array of semiconducting nanotubes (NTs) 
of length, 2h, placed in external electric field, F, parallel to their axes. By employing the conformal 
mapping, we find analytically the exact distribution of the induced charge density for arbitrary ratio 
of the gap, Eg, and the voltage drop, eFh, across the NT. We use the result obtained to (i) trace the 
narrowing of the neutral strip at center of the NT with increasing F; (ii) analyze the polarizability 
of the array as a function of E g ; (iii) study the field enhancement near the tips. 

PACS numbers: 73.63.Fg, 77.84.Lf, 79.60.Ht 
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I. INTRODUCTION 

During the past decade much progress has been 
achieved in growth of vertically aligned arrays of car- 
bon nanotubes (NTs). Vertical alignment is attractive 
for such applications as field emission displaysi^i^&i 
and chemical sensors^. Latest publications ^ 10 ' 11 : 12 re- 
port almost perfect alignment. Since recently, the efforts 
of researchers were also directed at synthesis of nanto- 
bube arrays aligned in a plane 13 : 14 i 15 : 16 i 17 ' 18 . The ob- 
jective of this effort is the synthesis of complex orga- 
nized NT structures for integrated molecular electron- 
ics devices and optoelectronic applications. Growth and 
patterning of in-plane oriented arrays require novel tech- 
nological approaches. A number of such approaches have 
been proposed in Refs. IT^IT5llT8l . and successfully imple- 
mented. 

A planar array of NTs is shown schematically in Fig. 1. 
In a typical experimental situation the distance, d, be- 
tween the neighboring NTs in the array is much smaller 
than the length, 2h, of the constituting NTs. This sug- 
gests that the response of the array to the external elec- 
tric field must exhibit a collective character. For verti- 
cally oriented arrays, the collective electrostatic response 
has been previously discussed in the literatur o 19 : 20 ' 21 . 
This response plays a crucial role in designing the NT- 
based field emitters. For planar arrays, being the emerg- 
ing area of research, their collective electrostatic proper- 
ties have never been addressed. 

Intuitively, it is obvious that responses of vertical and 
planar arrays to external electric field are dramatically 
different. Indeed, external field does not penetrate into 
a dense enough vertical array much deeper than the in- 
tertube distance 20 . In other words, to the first approx- 
imation, the tops of NTs constituting the vertical array 
can be considered as a metallic plate. By contrast, for a 
planar array, the force lines of external field can extend 
"above" and "below" the NT plane (see Fig. lb), thus 
causing the charge separation within each NT, consti- 
tuting the array. The collective character of the electro- 
static response manifests itself in the fact that the actual 
field, causing this separation in a given NT is strongly 
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FIG. 1: a) a planar array of semiconducting NTs in external 
field, F, is shown schematically; d is the distance between the 
neighboring NTs. As a result of charge separation, the posi- 
tively and negatively charged regions of length, h — ho, emerge 
near the tips. These regions are separated by neutral strips 
of width 2ho; b) Force lines of electric field above and below 
the array are curved as a result of the charge separation. In 
charged regions the force lines enter (exit) the NTs normally 
to the plane of the array. 



altered by the presence of the neighboring NTs. The 
resulting dipole moment, induced in an individual NT 
within the planar array, depends on the charge distribu- 
tion, the form of which is not obvious a priori. The situ- 
ation is even less trivial if a planar array is made of semi- 
conducting NTs. For a single semiconducting NT the 
charge separation happens for a field exceeding the crit- 
ical value Fth — E g /2eh. As F exceeds F t h, the charges 
occupy progressively larger portion of the NT; positive 
and negative charges are separated by a neutral strip 22 
of a length E g /eF. It is not obvious how, upon increasing 
F above F t h, the field-induced charge separation occurs 
in a planar array of semiconducting NTs. In particular, 
how the length of the neutral strip depends on the pa- 
rameters of the array, namely, the length and density of 
constituting NTs. 

The questions formulated above are addressed in the 
present paper. In Sect. II we extend the approach devel- 
oped in Ref. to the planar array of NTs. In Sect. Ill 
we demonstrate that for the planar array the Thomas- 
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Fermi equation can be solved exactly and find a simple 
analytical expression for the charge distribution along the 
NT for arbitrary external field. In particular, we find a 
universal relation between the length, 2h,Q, of the neutral 
strip and the applied field, F, and discuss qualitatively 
the limiting cases; Sect. IV. In Sect. V the exact form 
of the charge distribution is utilized to calculate analyti- 
cally various characteristics of the array, such as (i) elec- 
tric field inside the neutral strip, (ii) electric field outside 
the array, (Hi) enhancement of the electric field near the 
tips of NTs, and (Hit) polarizability of NTs in the array. 
Concluding remarks are presented in Sect. VI. 



II. THOMAS-FERMI EQUATION FOR THE 
PLANAR ARRAY 



Derivation of the Thomas-Fermi for the planar array 
is completely analogous to that for vertical array2°- (see 
also Ref. \2Si ) and goes as follows. The charge separa- 
tion in a given NT of the array occurs due to the ex- 
ternal field and the field, created by induced charges 
on all NTs. This separation proceeds until the local 
chemical potential assumes a constant value along each 
NT of the array. This local chemical potential of ra- 
th NT, p m (z), is related to the local charge density via 
the NT electronic spectru m 20 ' 22 as follows: p m {z) = 
sign(z)(E 2 /4 + [p m (z) / 'g} 2 ) 1/2 . Then the condition of 
the constant chemical potential can be presented in the 
form of generalized Thomas-Fermi equation, which for a 
given NT, m, in the array reads 
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eFz- / dz' y p n (z')<j)(z, z'\ fl TO - Rr, 
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(1) 



where \R n — R m \ is the distance between the centers of 
the NTs n and m. In Eq. ([1]) the dimensionless interac- 
tion parameter g is related to the electron velocity, Vq, 
in the NT as^ 2 - g = (2Ne 2 /Tre*hv ), where N is the 
number of conducting channels, which assumes the val- 
ues 2 or 4 depending on the spin and band degeneracy. 
This degeneracy is governed by the NT chirality. The 
potential, <fi(z, z'), in Eq. JT]) is given by 



(z,z';R) = 



V(z - z') 2 + R 2 ' 



(2) 



where e* = (e + l)/2 is the average dielectric constant 
between the substrate and the air. For a regular array 
we set all p n (z) equal to p(z) and introduce the NT den- 
sity, A/"o, as A/"o = 1/d. As the next step, we take the 
continuous limit of Eq. ([T]). In doing so, we first isolate 



the NT with n = and rewrite Eq. (JTJ) in the form 



eFz = sign(z) 



\ 



P(z) 



2C d p(z) 



+ ^ [ dz'p(z')S(z,z% 



(3) 



where R n — nd, and the kernel, S(z, z'), is defined as 



S(z,z') = £ 

njtO 



^{z-z>) 2 + R 2 n ^{z + z'f + Rl 



(4) 



Here we have used the fact that p(z) = —p(—z). The 
term 2Cdp{z) in the rhs of Eq. ([3]) comes from the term 
with n = in the sum |1| and describes the "self-action" 
of a given NT. For an isolated NT— this term has the 
form 2 hi(h/r)p(z), where r is the NT radius. In the pres- 
ence of neighboring NTs, the self-action is limited by the 
distance ~ d, so that Cd is equal to ln(d/r) analogously 
to the vertical array^S. Now taking the continuous limit 
in Eq. ([3]) reduces to replacement the sum over n in the 
kernel Eq. (j4|) by the integral, leading to the following 
integral equation 



eFz 



sign(z) 



\ 



2C d p(z) 



2Wn 



dz'p{z')\n{\z- z'\). (5) 



Rigorous justification of the continuous description is 
given in Appendix, where we show that the effect of dis- 
creteness of NTs in the array amounts to a small ~ d/h 
correction to p{z). 



III. SOLUTION OF THE INTEGRAL 
EQUATION 

A. Metallic NTs 



In the case of metallic NTs we set E g = in Eq. iJSJ . 
We will also neglect the self-action term, 2Cd, in the rhs; 
the validity of this step will be verified in the end of this 
Subsection. Then, for positive z, Eq. ([5]) takes the form 



eFz = ^ rdz'p(z') In 



(6) 



Note that Eq. ([6]) does not contain any small parameter. 
Indeed, upon introducing dimensionless charge density 
and coordinate as 



ee*F 



z 



(7) 
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Eq. assumes the following form 

z = I dz' p(z') In 
Jo 



takes the form 



(8) 



In Eq. ((8]) the relevant values of p, z', and z are all of 
the order of unity. It turns out, however, that, despite 
the absence of a small parameter, Eq. |(5J) can be solved 
analytically in the closed form. The solution reads 



7Tv/l-(f') 2 ' 



(9) 



In order to check that Eq. ([9]) is indeed the solution of 
the integral equation ijHJ), we substitute © into the rhs 
of Eq. ([8]) and perform integration by parts as follows 



dz' - = 

o ^1 - (i') 2 



In 



= _ 2; /W V- 1 "^! 2 =-22 



z — z' 

it/2 



(10) 



(5') 2 



(cos y>) 2 



(sin <p) 2 



where the integration in the rhs implies taking the prin- 
cipal value. The last identity in Eq. (| 10|) emerges upon 
substitution z — sin ip. Our key observation is that the 
integral in rhs is equal —n/2 for all \z\ < 1, so that p(z) 
in the form ([9]) is indeed the solution of Eq. (JSJ) . Return- 
ing to dimensional variables, we present the final result 
for the charge density distribution in the form 



p(z) 



ee*Fz 



nAfoVh 2 - z 2 



(11) 



Recall now, that this result was obtained neglecting the 
self-action. To verify the validity of this assumption, we 
estimate the ratio, 2Cdp(z)/ee* Fz, of the self-action term 
in Eq. ([5]) to the lhs. As follows from Eq. (TTTj) . this ra- 
tio is ~ Cd/Moh. Thus, the self-action can be neglected, 
if Afo > Cd/h, i.e., when the array is sufficiently dense. 
Although the self-action term is relatively small, it over- 
weighs the integral term near the tips, where the solution 
Eq. (|lip diverges. This suggests that the divergence in 
([TTjl must be cut off at (h — z) ~ h (£d/ 7lA/o) , where 
the the first term in the rhs of Eq. |(SJ| becomes compa- 
rable to the lhs. At the cutoff, the ratio (h — z)/h is 
~ (Cd/hhfo) <C 1, indicating that, in a dense array, the 
solution Eq. (JTTJ applies up to the close vicinity of the 
tip. 



B. Semiconducting NTs 

1. Reduction of the Thomas- Fermi equation to the 
Hilbert form 

Upon neglecting self-action, and in the limit of large g, 
the Thomas-Fermi equation for the charge distribution 



eFz - sign(z) -± = — / dz'p(z') ]n(\z - z'\). (12) 
2 e* J-h 



As was explained qualitatively in the Introduction [see 
also Ref. (|22T )]. in contrast to the metallic case, the solu- 
tion of Eq. ([12")) is singular. This is because the positive 
and negative charges near the tips, induced by the ex- 
ternal field, are separated by the neutral strip, caused 
by a finite bandgap, E g , with boundaries at z = ±/k). 
Within this strip we have p(z) = 0. As the external 
field increases, the strip narrows. Solving Eq. (fT2f im- 
plies finding both the form of p{z) outside the interval 
[—ho, ho] and the dependence of ho on the external field. 

Formally, with finite E gi the integral equation contains 
one dimensionless parameter, E g /eFh. In this subsection 
we demonstrate that the exact solution of Eq. (| 1 2[) can be 
obtained for arbitrary value of this parameter. As a first 
step it is convenient to differentiate both sides of Eq. (jX3J) 
with respect to z. This yields the following Hilbert-type 
integral equation 



2M 



(13) 



ft 



Our key observation is that the 2D geometry of the array 
allows one to employ the technique of conformal mapping 
in order to solve Eq. (fl"3|) . The choice of the appropriate 
conformal transformation is described below. 



2. Conformal Mapping 

The conformal mapping method can be applied for 
large g, when each of the regions —h < z < —ho and 
ho < z < h are equipotential. Then a general solution 
of the 2D Poisson equation with boundary conditions of 
equipotentiality on flat metal surfaces can be obtained 
from mapping of the boundaries of these surfaces on the 
real axis in the complex plane. What makes the case of 
semiconducting array nontrivial, is that the boundaries 
of equipotential regions have a doubly-connected geome- 
try. 

To gain an insight as to what is the appropriate form 
of the transformation for the geometry under study, it 
is instructive to start with an auxiliary problem of two 
isolated metallic plates in external field, occupying the 
intervals [—h,—ho] and [ho,h\. Two metallic plates cor- 
respond to the two cuts on the complex plane, as shown 
in Fig. 2. The general solution of Eq. (fT3|) for this ge- 
ometry can be found in the textbook o 24 ^ 25 ' 26 ' 27 . It can 
be expressed through the integrals along the closed con- 
tours L_ and L + , encompassing the cuts in the clockwise 
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Im(t) 



Re(t) 



FIG. 2: Contours of integration in Eq. (fT4)l , L- 
encompassing the cuts in the clockwise direction. 



and Lj 



The relation between ho and E g can be now established 
by substitution of the solution Eq. (fTB"|) with unknown 
constants A and ho into Eq. (fT3|) and equating linear in 
z and constant terms to eFz and E g /2, respectively. This 
yields A = —e*F/-K 2 N and the following transcendental 
equation for ho 



2eFh 



= E 



"[) 
/l. 2 



-h-flK 



/) 2 
h 2 



(17) 



where K(x) and E(x) are the elliptic integrals of the first 
and the second kind, respectively. In conclusion of this 



direction, in the following way 
A 



G=h /h 



y/(z*-h*)(h*-Z*) 




Vh 2 - t 2 y/t 2 - h 2 



z - t 



(14) 



3. Solution for the charge density distribution 

Although Eq. (fl4|) satisfies Eq. (fl"3|) . it cannot be ap- 
plied to the planar array of semiconducting NTs. The 
reason lies in the fundamental difference between the iso- 
lated metallic plates and array of semiconducting NTs. 
In the former case each plate respects electroneutrality, 

i.e., J_ h dzp(z) = 0, j h dzp(z) = 0, while the array is 
neutral only as a whole. Still Eq. (fT4)) can be modified to 
account for this difference. The modification amounts to 
adding a constant to the integral in Eq. IT4]) after which 
acquires the following form 



1 



^(z 2 - h 2 )(h 2 - z 2 ) 

fh 



xsign(^2A, ^VpS + C ; 



ho 



t 2 



(15) 



It can be checked by a direct substitution that Eq. (fT5|) 
is the solution of Eq. (fT5)) for arbitrary C. On the other 
hand, with nonzero C, the net charge in each of the 
intervals [— h, — ho] and [/io 5 ^] is finite. On the phys- 
ical grounds, the value of the constant C should be 
now determined from the condition p(ho) = 0, which 
fixes the ratio of the constants C and A at the value 
C/A = — (Tr/2)(h 2 — /i§). The meaning of this condition 
is that the charge density vanishes at the boundaries of 
the neutral strip. It turns out that the result Eq. (fTS)) 
can be greatly simplified. Namely, the integral can be 
evaluated analytically leading to the following simple ex- 
pression for the induced charge density 



p(z) = sign(z) 



ee*F / z 2 — h 2 . 



nAfo V h 2 - 



(16) 
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FIG. 3: (Color online) Curve 1: Dimensionless length of a 
neutral strip G(x) — ho/h is plotted versus dimensionless ra- 
tio x = Fth/F from Eq. (|17[) ; Curves 2 and 3 are the asymtotes 
plotted from Eq. (TK))) . 

Section, we note that conformal mapping was previously 
employed to find the positions of boundaries of the two- 
dimensional electron gas in a smooth external potential 28 
as well as to find the positions of edges of incompressible 
strips in the quantum Hall sample^. 



IV. DISCUSSION 

Eq. (fTT)) defines a universal dependence of the length, 
2ho, of the neutral strip on the applied field, F. This 
dependence can be parameterized as 



ho 



ko(!£), 



(18) 



where Fth = E g /2eh is the threshold value of F that 
causes the charge separation in the array. The dimen- 
sionless function, G(x) is shown in Fig. 2. It has the 
following asymptotes 



G{x) 
G(x) 



2^ 



x < 1; 



2 (1 - x) 
In (1 - x) 



{1-x) < 1. 



(19) 
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As seen from Fig. 2, the above asymptotes approximate 
G(x) very closely. The dependence G(x) should be con- 
trasted to the corresponding dependence for an isolated 
NT, where it has the formS 2 - Gq(x) = x. We note, that 
G{x) > Gq(x) for all x. In other words, the neutral strip 
for NT array is wider than that for a single NT. The dif- 
ference is especially pronounced for strong fields, corre- 
sponding to x <c 1, where we have G(x)/Gq(x) ~ x~ x / 2 . 
The origin of such a difference lies in the fact that the 
array screens the external field much more efficiently 
than a single NT. As a result, for strong fields, when 
h <C h, the field inside the strip is ~ (h /h) times 
smaller than F. Then the condition that the poten- 
tial drop within the neutral strip is equal to E g can be 
presented as eho(hoF/h) ~ E g . This condition yields 

1/2 

ho ~ (Egh/eF) , which is in agreement with the small- 
x asymptote G(x) ~ a; 1 / 2 of Eq. (TlT)|) . 

For an individual NT, the induced charge density 
changes linearly near the edges of the neutral strip 22 . 
By contrast, for a 2D array, p(z) described by Eq. (fT6|) 
exhibits a singular behavior, p{z) oc ±y/ z ± ho, near the 
boundaries of the neutral strip. 

The principal assumption that allowed us to find the 
charge distribution in semiconducting NT analytically is 
that in Eq. (12]) we had replaced y/(E g /2) 2 + (p/g) 2 in 
Eq. ([5]) by E g /2. Using the result Eq. (fTB|) . we can now 
check the validity of this assumption. It is easy to see that 
the replacement is justified when the external field is not 
too strong, so that the dimensionless ratio x — F t h/F in 
Eq. (TT5)) must exceed the value (gNoh) -1 <C 1. 



IMPLICATIONS 



Exact results Eqs. (fTTj) , (TTS)) for the charge distri- 
bution in metallic and semiconducting arrays, obtained 
above, enable us to calculate analytically various observ- 
able characteristics of the array. In this Section we will 
focus on the following characteristics: electric field dis- 
tribution inside and outside of the array, and collective 
polarizability of the array. Concerning the electric field 
distribution, the questions most interesting for applica- 
tions are the distribution of electric field within the neu- 
tral strips and behavior of electric field near the NT tips. 
The fist question is important for electrooptical measure- 
ments, while the second question is crucial for the field 
emission. In addition, it is instructive to study quantita- 
tively how the charges, induced in the NTs by external 
field, alter this field in the regions above and below the 
array. In particular, the interesting question is how the 
spatial decay of the electric field component normal to 
the array depends on the position, z, along the NTs. 
These questions are addressed in subsections that follow. 



A. Electric field inside the neutral strip 

Electrostatic potential inside the neutral strip is given 
by superposition of the external potential, eFz, and the 
induced potential 



(j) md {z) 



dz' p(z') In 



1 1 o 



(20) 



where z < ho- Upon substituting the exact form of 
the charge density, Eq. (T5J) , in to Eq. PO)) the integra- 
tion can be performed analytically. We will present the 
result for the tangential component of the net electric 
field which is the sum of the external, F, and induced, 
pmd _ Q(j) md ( z y Q z , fields inside the neutral strip, 



f; 



F. 



ind 



F = F\ 



[ho(F)]' 



h 2 



(21) 



where the dependence ho(F) is defined by Eq. §17\ and 
plotted in Fig. 3. In the limit of weak external field, close 
to the threshold value, F t h, Eq. pTjl yields a natural re- 
sult that Ff w F. This is because the induced charge 
occupies only a small portion, (h — ho) <C ho, of the 
NTs. In the opposite limit of the strong external field, 
F F t h, the field, Ff, at the center of the NTs is equal 
to Fho[F]/h <C F, i.e., it is strongly suppressed, and 
falls off towards the ends of the strip. Note, that in our 
qualitative derivation of the width of the neutral strip 
we have used the estimate Ff ~ Fho/h for the field in- 
side the strip. Therefore, Eq. (|21| justifies this estimate. 
The reason why we focused on the tangential compo- 
nent, F**, is its possible relevance for the experiment^ on 
electroabsorption of light. It might seem that, with pho- 
ton energy ~ leV much bigger than eFh, the large-scale 
nanotube geometry is not important. This, however, is 
not the case. The reason is that the dipole moment of 
the many-body optical transition is directed along the 
tube^ii^ 2 .^^^. In analysis of the data in Ref. it 
was assumed that in the change, [a(u,F z ) — a(u>,0)], of 
the absorption coefficient, F z is the external field. Mean- 
while, the optical transition is modified by the acting 
field, Ff. Then Eq. J2T]) suggests that the acting field 
can be much smaller than the external field and, more- 
over, behaves as Ff cx F 1 / 2 due to the collective screen- 
ing. 



B. Electric field outside the array 

Consider a point with coordinate, z, located at a dis- 
tance, H, above the array. Similarly to Eq. (|2"Tj) , the 
parallel and perpendicular components of the acting field 
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can be evaluated from 



F z (z,H)=F- — / dz'l dx- 



/• ;•-•//; 7777" I " !; " 



(z' - z) 2 + x 2 + H 2 
(z' - z) 2 + x 2 + H< 



We will restrict ourselves to metallic NTs. In this case 
the integration in Eq. (|22|) can be performed analytically, 
yielding 




F z = 



/2A 
V2F 
~7I 



(z 2 + H 2 )VX + A-h 2 (h 2 + H 2 - z 2 ), 



\] (z 2 + H 2 )VI - A + h 2 (h 2 + H 2 - z 2 ), 



(23) 



where we had introduced an auxiliary function A(z, H) 
defined as follows 



A(z, H) = (H 2 + z 2 - h 2 



Ah H 



(24) 



Eq. ([2"3"| allows one to trace how the tangent component, 
F z , grows from zero to F with increasing H, while the 
normal component falls off from 2nJ\fQp(z) to zero with 
increasing H . The characteristic scale of change for both 
components is H ~ h. This behavior is illustrated in 
Fig. 4. 



H/h 



0.5 



1.5 



FIG. 5: (Color online) Dimensionless tangent component, 
F z /F, of the net electric field is plotted from Eq. ([23} vs. 
dimensionless distance, H/h, from the array for various val- 
ues of coordinate, z, along the array. The plots are presented 
from bottom up: z/h = 0.1; 0.7; 0.9; 0.95. 



and z/h — 0.9 intersect at H/h « 0.3677 (see inset 
in Fig. 5). This implies that the derivative dF±/dz 
is zero at this point. On the other hand the relation 
dF±/dz = dF z /dH suggests that F Z (H) must have a 
maximum at H/h w 0.3. As seen from Fig. 5, this is 
indeed the case. It is also seen from Fig. 4 that, for z/h 
approaching 1, the field F±_ grows in the vicinity of the 
plane of the array H — > 0. This is the manifestation of 
the enhancement of electric field near the tip. We address 
this issue in more detail in the next subsection. 




_| — , — , — , — , — , — , — , — , — , — , — . — H/h 

0.5 1 1.5 2 

FIG. 4: (Color online) Dimensionless normal component, 
F±/F, of the net electric field is plotted from Eq. (|23p vs. 
dimensionless distance, H/h, from the array for various values 
of coordinate, z, along the array. The plots are presented from 
bottom up: z/h = 0.3; 0.6; 0.85; 0.9. Two insets illustrate 
that the curves F±(H) for z/h — 0.8 and z/h = 0.9 intersect 
each other. The actual point of intersection (not showed) is 
at H/h « 0.3677. 

We see that, as could be expected, F± falls off mono- 
tonically with increasing the distance, H, from the ar- 
ray. Less trivial is that the curves for z/h = 0.85 



C. Field enhancement near the tips 

Consider a point, zq, on the axis of one of the NTs 
outside the array, zq > h. The field in this point has 
only z-component given by 

FT d -^i h d ZP {z)(^- —\ (25) 

ee* J \ z + z z -zj 

With metallic charge density distribution Eq. (fTTj) the 
integration in Eq. ([23)) can be performed analytically. We 
will present the result for the field enhancement factor, 
(3 = F" Ld /F, which is the conventional characteristics of 
the field near the tip 



We can now compare the dependence (3{zq) for the pla- 
nar array and for an isolated NT. For the isolated NT 
we have^ 2 - ft(za) cx h/(zQ — h) (in the present nota- 
tions). By contrast, Eq. ([2"6l) yields a slower decay, 
/3(zq) ~ h 1 / 2 /(z — h) 1 ! 2 . On the other hand, the mag- 
nitude of the enhancement factor is much weaker for the 
planar array than for the isolated NT. Obviously, this 
reduction is due to the collective screening. 
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D. Induced dipole moment 

The dipole moment, P, of a given NT within the array 
is defined as 

P(F) = 2 f dz zp(z). (27) 

Analytical evaluation of the integral in Eq. (|27|) . with 
the charge density fH]), yields the following expression 
for the dipole moment 

P(F) = ^{h 2 -{h (F)} 2 ). (28) 

It is instructive to consider the limiting cases of weak 
and strong fields. Using the low-field and strong-field 
asymptotes of h (F) given by Eq. (fl9|) . we obtain 

Note that in the limit of low fields, (F — F t h) <C F t h, the 
behavior of P{F) is singular, unlike the individual NT, 
where one has^ P(F) oc [1 - F/F th ] . This singularity 
is the manifestation of the singular behavior of the charge 
density, p(z) oc \J z — ho, at the boundary of the neutral 
strip. By contrast, for the individual NT, this behavior 
is regular, namely p(z) oc (z — ho). 

In the strong-field limit, the contribution of the neutral 
strip to P(F) constitutes a relatively small correction 
~ Fth/F. Essentially, the polarizability x = P{F)/F of 
a given NT is ~ eh 2 /Mo, as seen from Eq. |[2"B]). On the 
other hand, the value of x f° r the individual NT can be 
estimated from the textbook problem 3 ^ of polarization 
of an ellipsoid with large ratio of axes in electric field 
directed along the major axis. With logarithmic accuracy 
one has, x ~ e ^- 3 - Thus, compared to the individual 
NT, the polarizability of a given NT within the array is 
suppressed by a large factor Af h. This conclusion can 
be reformulated as follows. An element of the array of 
a width h contains A/q/i NTs. The integral polarizability 
of this element is ~ (Noh)(eh 2 / Wo) = eh 3 , i.e., it is of 
the same order as polarizability of the individual NT. 
Therefore, the magnitude of suppression of x within the 
array, can be understood from the fact that, within the 
array, the constituting NTs "talk to each other" only if 
their separation is smaller than h. 

VI. CONCLUDING REMARKS 

Throughout the paper we assumed that the NTs in 
the planar array are parallel to each other. Note, that 
recently a different class of planar NT arrays in the form 



of a random network, residing on the polymeric sub- 
strate, and connecting the source and drain electrodes 
has been discussed and realized experimentall y 37 ' 38 ' 39 . It 
was demonstrated that such a network can emulate the 
channel of the field-effect transistor. Since in our con- 
sideration of the planar array was carried in continuous 
limit, the results obtained are relevant to the random 
planar dense arrays. 

Our main results are Eqs. (fTTj) . (fT6|) describing the 
distribution of induced charge in metallic and semicon- 
ducting planar arrays, respectively. These results were 
obtained in the continuous limit, i.e., we replaced the ar- 
ray of individual NTs with a two-dimensional continuous 
layer of a zero width. Such a replacement is justified with 
high accuracy as long as the NT length, h, exceeds the 
separation, d, between the neighboring NTs. More quan- 
titatively, as shown in Appendix, the relative correction 
to the induced charge distribution is ~ d/h. Similarly, 
the randomness in distances between the neighboring 
NTs is not important for the planar array. This should 
be contrasted to the infinite vertical array 2 ^, where the 
distance between neighboring NTs governs the penetra- 
tion depth of the external field. By contrast, the con- 
tinuous description adopted in the present paper applies 
with high accuracy, since, in the planar array, the in- 
duced charge extends over the entire length, 2h, of the 
NTs. This length the is much bigger than Af ~~ ■ The dis- 
creteness of the array manifests itself only in the close 
vicinity of the tips, where it cuts off the singularity 
p(z) oc (h — z)~ 1/2 in the distribution of the induced 
charge. 

It is instructive to trace how the results obtained in the 
present paper for the planar array evolve when one con- 
siders the vertical array of a finite thickness, W. Consider 
first the most interesting case of intermediate thickness, 
d <C W <SC h. The modification of our result Eq. (TT]) to 
this case can be understood from the following reason- 
ing. Each square with the side W contains (W/d) 2 NTs. 
This square can be viewed as a single "combined" NT. 
Then the charge density in this combined NT is given by 
Eq. (fTTj) with J\fo replaced by W^ 1 . To find the charge 
density per a single NT, this result should be divided by 
(W/d) 2 . Thus we conclude, that modification of Eq. (fTTj) 
to finite W < h amounts to the replacement of A/q -1 by 
d 2 /W, so that the induced charge density falls off as 1/W 
with increasing the width. However, the applicability of 
Eq. dTTJ) to the finite-width array is limited by the size 
of a "combined" NT, i.e., (h — z) must exceed W. As 
W exceeds the NT length, this condition cannot be met, 
and the charge distribution changes dramatically from 
the power-law behavior to the exponential screening 20 . 

Note that our basic integral equation, Eq. ([5]), express- 
ing the fact that electrochemical potential is constant 
along each NT of the array allow the exact solution us- 
ing the technique of conformal mapping only when the 
"self-action" [second term in the rhs of Eq. |5])] is ne- 
glected. This step is well justified for a dense array. In 
the electrostatics of inhomogeneous fractional quantum 
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Hall liquid a i the charge distribution is described by 
the equations similar to Eq. ([5]), i.e., with logarithmic 
kernel, but with self-action playing an important role. 
Then these equations can be solved only numerically. 



Therefore, after performing the summation over n in 
Eq. J4} and using Eq. ([35]). the rhs of Eq. © will ac- 
quire the form 
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VII. APPENDIX 

Below we present a rigorous derivation of Eq. ([5]) for a 
periodic array of NTs. We start from the identity 



1 



\Jx 2 + b 2 



-J dq K (\x\ q^j exp(iqb) . 



(30) 



where Kq(x) is the modified Bessel function of the second 
kind. To use this identity in Eq. ([3]), we set b — > R n = nd. 
Then the summation in the integrand in the rhs can be 
readily performed; it transforms the sum over n in Eq. Q 
into the following sum of the <5-functions 



n p \ / 



(31) 



where p assumes all integer values. As the next step we 
isolate the term with p — in the kernel Eq. ^ from all 
other terms with nonzero p. Then we make use of the 
fact, that 

K Q (x) =-7 + ln2-ln|x|+0(a; 2 ), x < 1, (32) 

where 7 is the Euler constant. Then, for the term with 
p = 0, we get 



- d lim [K (\z - z'\q) - K (\z + z'\q)] = ~ In 



S(z,z') = -In 



(35) 



First term in 



2np\ 



d ) U V~ ' 1 d J 

35]) describes the continuous limit; it re- 
produces the second term in the rhs of Eq. ([5]) and comes 
from p = in Eq. (|31[) . The remaining sum over p gives 
rise to the "self-action" term in Eq. ([5]) . The easiest way 
to see this is to replace the summation over p by integra- 
tion dp, which would immediately yield 



S (z,z') = <t>{z-z')-4>{z + z') 



(36) 



for the second part of the kernel l|35[) . Note that Sq(z, z') 
is nothing but the kernel for a single N T 20 ' 22 . However, 
the replacement of the ^ by J dp is justified only when a 

large number of terms contribute to the sum. This is the 
case only when the condition \z — z'\ < d is met. Then, 
upon substituting Q36[) into (|31|) and restricting integra- 
tion over z' to \z'—z\ < d, the self-action term in Eq. ([5]) is 
recovered. Concerning the second part of the kernel 



> d, it is dominated by the terms with p = ±1, 
1/2 



at \z' 

which behave as (d/\z~z"\j exp j— 2n\z— z'\/d\. Con- 
tribution coming from this part to the integral Eq. (|34| 
is ~ p(z), which is parametrically smaller than the self- 



(33) action term 2Cdp(z). 
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